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Steinhausen ' s  hypothesis,  that  a  simple  torsion-pendulum  model  adequately  describes  the  dynamics  of 
semicircular-canal  action,  has  received  its  sharp  of  criticism-  This  is  primarily  because  of  the  fact 
that  the  time  courses  of*  afferent  discharges  and  vestibular  reactions  to  angular-accel erati on  stimuli 
only  approximately  relate  to  the  time  course  obtained  as  a  solution  to  the  trrsion-pendul um  model  for  the 
same  angular  accelerations.  Thus,  there  is  a  need  t  o  determine  a  generalized  transfer  function  whi^h 
delineates  those  components  that  provide  meaning  ti,  Uteinhauspn 's  hypothesis  from  those  which  describe 
other  phenomena  such  as  "adaptation". 


INTRODUCTION 


Rotational  head  movements  stimulate  the  sensory  end  organs  of  the  semicircular  finals  of  the  inner  ear, 
causing  compensatory  eye  movements  (nyst agmus ) ,  appropriate  to  the  plane  of  rotation,  via  a  vestibtilo- 
ocular  reflex  arc.  However,  the  sensory  epithelia  ..f  the  semicircular  canals,  the  erlstae  ampul  lares, 
are  relatively  inaccessible  for  direct  functional  .study  during  normal  physiological  head  movements. 

In  *he  first  place,  they  lie  deep  within  the  skull,  enclosed  by  spongy  and  petrous  bone  and  the  osseous 
labyrinth.  Secondly,  any  technique  which  exposes  the  hair  cells  or  their  v>erv«*  supply  for  neurotiological 
experimentation  necessarily  disrupts  the  vascular  supply  to  the  neuroepithelium  and  interferes  with  the 
ionic  balance  between  the  perilymph  and  endolymph.  Thirdly,  even  when  the  membranous  ampullae  are 
exposed,  it  is  difficult  to  perceive  the  extent  to  which  the  normally  invisible  cupula  covers  the 
crista  and  fills  the  ampulla  of  the  duct.  (It  is  the  cupula  which  couples  the  mechanical  energy  of  head 
motion  to  the  receptor  hair  cells  of  the  crista  for  transduction  into  neural  impulses  at-  the  afferent 
terminals.)  Thus,  it  is  not  surprising  that,  for  over  a  century,  the  dynamic  response  characteristics 
of  the  semicircular  canals  have  not  been  completely  elucidated  despite  continuous  investigation  by 
many  vestibular  physiologists. 

Historically,  until  the  time  of  Fluorens  (lA),  the  semicircular  canals  were  associated  with  the 
perception  of  sound.  By  sectioning  the  membranous  semici reti  1 ar  canals  in  the  pigeon  and  the  rabbit, 

Fluorens  established  that  eye,  head,  and  body  movements  always  occur  in  the  same  plane  as  that  of  the 
injured  canal,  floltz  (l6),  also  investigating  the  pigeon,  first  associated  these  canals  with  bodily 
equilibrium.  Breuer  (5),  Crum  Brown  (9),  and  Mach  (.’6)  suggested  that  the  crist.ae  ampullares 
respond  to  rotatory  accelerations  of  the  head.  More.,ver,  they  also  hypothesized  that  the  motion  of 
the  endolymph  within  the  canal  was  re-sponnible  for  eliciting  the  neural  activity  in  the  nerve  endings 
(hydrodynamic  theory).  (Mach,  however,  later  appears  to  have  favoured  a  hydrostatic  theory,  i.e., 
that  pressure  is  the  adequate  stimulus — see  refs.  7f*,  Pfi.  )  Evald  (li),  using  an  "hydraulic  hammer", 
produced  artificial  endolymphatic  movements  in  the  canals  of  pigeons  and  demonstrated  that  the 
direction  of  head  nystagmus  was  dependent  upon  the  direction  of  endolymph  flow.  Evald's  "law"  was 
placed  on  a  firmer  footing  when  Lovenstein  and  Sand  (?5)  established  that  the  action-potential 
(impulse)  frequency  in  fibers  innervating  the  cristae  varied  according  to  the  direction  of  angular 
movement  of  the  head. 

In  order  to  describe  the  dynamics  of  the  semicire-ijar-canal  system,  rte inhausen  (76,  ^7),  and, 
subsequently,  F^mond ,  Oroen  and  .Tonkees  (10),  proposed  the  so-called  torsion-pendul  um  model.  In 
this  model,  the  Instantaneous  angular  deflection  of  the  cvjpula,  E(t),  ennped  by  an  instantaneous 
angular  acceleration  of  the  head,  o(t),  is  governed  by  the  linear  differential  equation 


0  d  +  n  -£-(^1  ♦  At f t  )  =  YOa(t), 

dt?  dt 


(1) 


where  0  is  the  effective  moment  of  Inertia  of  the  endolymph  in  the  semicircular  canal ;  H  is  the 
viscous  damping  momeni  per  unit  angular  velocity  of  endolymph  relative  to  the  wall  >T  the  membranous 
canal;  A  is  the  elastic-restoring  moment  per  unit  angular  displacement  of  the  end-lymph  relative  to 
the  membranous  ampul lary  walls;  and  t  is  a  constant  of  proportionality  relating  ratio  volumetric 
displacements  (between  the  cupula  and  the  canal  endolymph).  The  transfer  function  pq,  (l)  can 
be  written  as : 
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;ne  way  t,.  1  ■  the  adequacy  •?'  this  model  is  t.>  r*»c  the  neural  activity  <f  the  primary  afferent 
fibers  innervating  t  he  semicircular  canals  a.;  they  resjiond  t.<>  appropriate  stimulation.  This  has  been 
done  f<-r  a  variety  -*f  specie;:  ,f  animals  (2,  1,  12,  17,  IP,  *!3,  2u ,  33,  3l< ,  ■!>,  3ft)  and  the  model  has  been 
found  to  he  deficient  .,  The  present  set  of  experiments  was  designed  t'_  study  the  mechanoneural  response 
ohuruoterist lea  of  primary  afferent,  semicircular-canal  units  in  the  pigeon,  in  order  to  determine 
whether  r  not  a  generalized  transfer  function  e,.uid  be  obtained  which  would  also  describe  similar 
neuro dynamic;:  in  other  species.  (In  testing  to  this  model,  the  assumption  is  made  that  the  afferent 
activity  is  proportional  to  the  cupular  displacement.) 


METHODS 

White  King  pigeons  ( Ool umba  H  via )  were  surgically  prepared  for  microelectrode  recording  from 
peripheral  units  which  innervated  the  semicircular  canals  (primarily,  the  anterior  semicircular  canal). 
The  animal  (with  head  immobilized)  was  oriented  on  board  a  rotatory  device  so  that  the  center  of  its 
head  plane  was  coincident  with  t.he  plane  of  rotation  (according  to  the  type  of  canal  being  recorded 
from ) . 

The  main  rotatory  sequence  consisted  of  a  scries  of  sinusoidal  angular  accelerations, 

aft)  =  a  sin  Pirft,  (3) 

m 

which  were  delivered  to  anesthetized  preparations  at  frequencies,  f,  from  0.01  to  10  Hz  with  peak 
angular  accelerations,  am,  of  2.0,  1».0,  8.0,  12.0,  and  20,0  degrees /s^  (t  ^  real  time  in  Eq.  (3)). 

The  single  unit  neural  activity  was  amplified,  displayed  on  an  oscilloscope  and  recorded  on  magnetic 
tape  according  to  conventional  techniques.  (One  channel  of  the  tape  recorder  was  used  for  voice 
commentary;  another  for  the  stimulus  reference  signal.) 

The  stimulus  reference  signal  was  used  to  trigger  a  physiological  signal  analyzer  (Nicolet 
Instruments  Inc.)  t.o  count  and  store  the  taped,  entrained  action  potentials  in  preselected  appropriate 
time  periods  (bins)  for  up  to  U096  sequential  bins  (see  Fig.  1  for  typical  binned  response). 

Fourier  techniques  were  used  on  the  binned  neural  data  to  obtain  the  magnitude  of  the  peak  amplitude  of 
the  response  and  the  temporal  (phase)  relationship  between  the  angular  acceleratory  stimulus  and  the 
fundamental  component  of  the  neural  response. 

Amplitude-  and  phase-values  were  usea  in  a  curve-fitting  program  on  a  PDF'-ll/bO  minicomputer 
(Digital  Equipment  Corp. )  to  provide  a  distinct  mathematical  expression  for  the  best-fitting  transfer 
function  for  a  linear  system.  (The  steady-state  frequency  response  Tor  a  linear  system  to  a  sinusoidal 
input  angular  acceleration  may  be  found  from  gain, 

!0(f)|  =  ((He  0(f))L>  +  (Im  (Hf))'-’}5*,  (M 

and  phase , 

0(f)  =  tan-1  dm  <i(f)/Re  G(r)),  (5) 

spectra  ( Bode  plots — see  ref.  30),  where  Re  0(f)  and  Im  0(f)  are  the  real  and  imaginary  parts  of  G(f), 
respectively. ) 

The  method  employs  techniques  in  non-linear,  least  squares  approx. mations  and  is  applicable  to 
both  high-  and  low-order  transfer  functions  0»).  Ori  the  final  Interpol*  tlon  in  the  program,  the  least 
squares  error,  LSE,  of  the  best-fitting  transfer  function  is  obtained  as 

L&E  *  YTY,  (<>) 

where  Y  is  a  residuals  matrix  (Y*^=  transpose  of  Y)  which  is  made  up  of  error  terms  that  express  the 
differences  between  the  experimental  data  and  the  model  which  is  to  be  fitted.  The  mean  square  error, 

M15E  *  LSE/(2L  -  S),  CO 

compares  the  goodness-of-fit  of  the  derived  transfer  function  to  that  of  the  experimental  data  (2L  £ 
number  of  real  and  imaginary  components  of  the  data  points,  and  i  number  of  parameters  in  the  transfer 
function ) . 

A  more  detailed  description  of  the  methodology  may  be  found  in  a  special  monograph  (7)  anu  elsewhere 
(8,  21). 


result:4  and  discussion 

The  simplest  transfer  function  that  fitted  the  data  for  all  units  ia  of  the  form. 
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(0) 


where  s  i «  Cract  i 1  —  -rder  differential  operator  with  0  <  k  <  1,  and  n  is  a  gain  constant  with 
units  in  impulu*-:--s~1/ii»  grct-:-,  (.’1).  The  Bode  pints  for  four  of  the  units  are  shown  in  Figs.  2-5, 

together  with  their  Best -fitting  ami  the  best-fitting  torsion-pendulum  model 


(In  rigs.  2-5,  the  effects  of  the  To  —  2.0  ms  (vide  infra)  contributes  no  more  than  ll  to  the  gain  spectra 
between  f  -  0.01  and  10  Hz. ;  consequently,  the  single-pole  transfer  function,  0"(s),  was  an  adequate 
representation  of  the  torsion-pendulum  model.)  As  is  evident  i»  the  plots  in  Figs.  2-5,  G'(s)  is  a 
much  be*  ter  *'it  to  the  data  than  is  G"(s)  (of.  the  M3F.  values  for  the  two  models);  and,  in  particular, 
the  fit  appears  to  improve  with  increasing  k. 

K 

Wnat,  then,  is  the  significance  of  k,  or  better  still,  s  ?  Elsewhere,  it  is  shown  that 


M 

*  K  n 

i  =  l 


(10) 


where  K  and  Bj  are  constants,  H  s  *  (in  theory,  but  finite  when  fitting  Eq.  (10)  via  a  digital 
computer),  and  and  tj  j  { Cq  :  0)  are  time  constants  (21),  Interestingly,  when  M  =  1,  Eq.  (10) 
becomes 


(11) 


which  has  previously  been  defined  as  the  transfer  function  of  the  adaptation  operator  (27,  1*0).  Thus, 
appears  to  be  a  form  of  adaptation.  Work  by  Thorson  arid  Riederman-Thorson  (39)  suggests  that 
sk  represents  a  distributed  relaxation  process  which  is  inherent  in  the  sensory-adaptation  mechanics 
of  Li  mu  1  us  photoreceptors,  vertebrate  retinal  receptors,  chemoreceptors ,  and  other  mechanoreceptors . 
Investigations  by  Taglietti,  Rossi  and  Casella  (30)  further  suggest  that  s*1  likely  represents  a 
relaxation  phenomenon  consisting  of  a  time-varying  intracellular  electrogenic  process,  the  com¬ 
ponents  of  which  are  summed  with  the  generator  potential  in  the  receptor  hair  cell. 


The  coefficient,  of  variation,  CV ,  was  determined  as  the  ratio  of  the  standard  deviation  of  intervals 
to  the  mean  interval,  as  obtained  from  int.erspi ke-interval  distributions  of  spontaneous  single  unit 
activity.  When  a  regression  of  CV  on  k  was  made  for  20  units,  a  statistically-signif leant  product- 
moment  correlation  (r  =  0.301*,  F’  <  0.05)  was  obtained  (21).  Thus,  the  larger  the  CV  is,  the  larger  the 
value  of  k  and,  consequently,  the  amount  of  adaptation.  Other  work  by  Goldberg  and  Fernandez  (15).  in 
squirrel  monkeys,  shows  that  the  CV  Is  statistically  correlated  with  semicircular-canal  afferent 
fiber  conduction  rates.  The  thicker  fibers  have  faster  conduction  rates  and  larger  CVs.  Together, 
these  findings  suggest  that  sensor 1  adaptation  phenomena  are  directly  dependent  on  the  innervation 
pattern  of  the  afferent  fibers. 

The  transfer  function,  G*(s),  differs  also  from  that  of  Ci"(s)  in  that  is  not  single-valued  as 
it  is  in  the  torsion-pendul um  model;  rather,  it  is  unit  dependent,  taking  on  values  from  t^  *  1*.1*5  to  22.17 
s  (mean  ±  GEM  *  10.2l»  ♦  1.20  s)  (21).  (in  fitting  Eq.  (0),  the  coupling  between  k  and  would  account 
for  some  of  the  five-fold  range  of  values  that  were  determined  for  tL.  However,  there  is  sufficient 
indication  from  other  studies  (32)  that  the  response  dynamics  of  small  groups  of  contiguous  hair  cells 
are  quite  different  from  those  of  other  groups.)  Realizing  that  the  nair-cell  tufts  are  stiff  (lU), 
that  their  lengths  vary  according  to  their  position  on  the  sensory  epithelium  (crista)  (22),  that  the 
number  and  thickness  of  the  stereocilia  can  be  variable  (22),  and  that  the  mechanical  properties  of  the 
cupula  are  not  necessarily  uniform  across  the  crista  (29),  then  it  is  plausible  that  could  have 
a  regional  distribution. 


The  form 

W  1  C«K 

(x^nHigS*1)  •  (l2) 

or  sometimes, 

G(s)(t  s+1 )  (13) 

m 

was  fitted  to  published  afferent-response  data  in  the  squirrel  monkey  (Salmiri  sciureus)  (12),  the 
frog  (Rana  esculents  and  R.  temporarla)( 3) ,  the  gerbil  (Meriones  ungulculatusT  ( 35 ) ,  and  the  guit&rfish 
(Rhlnobatos  productus)  (33).  The  parameter  xm  =  high  frequency  time  constant  which  results  from  both 
the  displacement  and  the  rate  of  displacement  of  the  cupula.  Such  a  term  has  been  obtained  from 
analysis  of  vest i bul ar-dri ven  eye  movements  In  man  (Benson  and  Sternfeld,  cited  in  ref.  l),  and 
primary  afferent  canal  responses  in  the  squirrel  monkey  (12)  and  the  elasmobranch  fish  (2U).  Table  1 
lists  k,  T.,  t ,  and  for  these  four  species  and  for  unit  responses  to  white-noise  stimuli  obtained 
from  pigeons  that  were  primarily  enc^phale  isol£  preparations.  For  all  species  listed  in  Table  1,  the 
MGE  using  Eqs.  (12)  or  (13)  was  comparable  to  or  smaller  than  that  obtained  with  other  models.  Further 
details  are  given  elsewhere  (0). 


!  .lift. i  >ti  (Eqs.  M:'>  nr  (13))  describing  jemirircular-ranaJ 
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TABLE  1:  Parameters  >.  f  a  genera 
dynamics  in  five  self  1 


'  '•  Tm 
sP*cies  i  ms  1  (s) 

frequency  range, 
f  •  H?  ) 

Squirrel  monney 

(a)  "regular"  ur.i’. 

(b)  "irregular"  «.  •* 

Gerbi  I 

(a)  CV  <  0  1 

(b)  CV  >  0.  I 

Fmg 

Guitarfi  sh 


0.0125-R 

0.0125-8 


0.01-5 

0.01-5 


0.0125-0.5 

0  ,w2-U 


0.5-16 


In  Table  1,  the  frequency  ■••» 
adequate  value  of  Tg.  In  gene/  ■» 
from  the  biophysical  proper*  i«v 
of  the  membranous  vestibular  at  t  *» 
Money  and  colleagues  (U)  f  r  *;• 
(f  =  l/2tiTj,)  of  80  Hz  in  the  tv* 
from  canal'  afferents  in  the  g  .  : 
given  evidence  that  5.70  ^  t  • 
fitting  to  afferent  data  fn  t  k<«* 


*  »  restricted  in  the  frog  and  guitarfish  to  utilize  an 

•  •  (s  0/11)  in  the  fits  have  been  predetermined  indirectly 

. ymph  and  the  dimensions  of  the  pertinent  anatomical  features 

•  .* a . > Jt*  irj  =  2  ms,  which  was  determined  biophysical ly  by 
.  .cgf-it;  that  there  should  be  an  upper  break  t'requency 

r«'..pi  ii.-.e  dynamics  in  semicircular-canal  afferents.  Recording 
.  .t.  , .  n'k  frequency  up  to  70  Hz,  Hartmann  and  Klinke  ( l8 )  have 
•t  id  the  same  range  as  those  determined  empirically  for 

i.  ,  .’1 ,  <5). 
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